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Abstract. In this paper, we consider a competition model between n species 
in a chemostat including both monotone and non- monotone response functions, 
distinct removal rates and variable yields. We show that only the species 
with the lowest break-even concentration survives, provided that additional 
technical conditions on the growth functions and yields are satisfied. LaSalle's 
extension theorem of the Lyapunov stability theory is the main tool. 



1. Introduction 

In this paper wc study the global dynamics of the following model of the chemo- 
stat in which n populations of microorganisms compete for a single growth-limiting 
substrate: 

s'(t) = D[S°-s(t)]-YZ:iMS(t))xi(t) 

x'M = ]pi(S(t)) - Di]xi(tl i = l---n, 

where 5(0) > and £j(0) > 0, i = 1 • • • n and S°, D and Di are positive constants. 
In these equations, S(t) denotes the concentration of the substrate at time t; Xi(t) 
denotes the concentration of the ith population of microorganisms at time i; fi{S) 
represents the uptake rate of substrate of the ith population; pi (S) represents the 
per-capita growth rate of the ith population and so the function yi(S) defined by 

Vi(S) = m 

is the growth yield; S° and D denote, respectively, the concentration of substrate in 
the feed bottle and the flow rate of the chemostat; each Di represents the removal 
rate of the ith population. We make the following assumptions on the functions pi 
and 

• Pi, fi : R+ — > R+ are continuous, 

• P»(0) = /i(0) = and for all S > 0, Pi (S) > and fi(S) > 0, 

For general background on model (|l.ip . in the constant yield case yi(S) = Y- L , 
the reader is referred to the monograph of Smith and Waltman |13j . Following 
Butler and Wolkowicz [3], we make the following assumptions on the form of the 
response functions pc there exist positive extended real numbers Xi and fj,i with 



Date: February 25, 2010. 

2000 Mathematics Subject Classification. 92A15, 92A17, 34C15, 34C35. 

Key words and phrases, chemostat, competitive exclusion principle, Lyapunov function, global 
asymptotic stability, variable yield model. 



2 



TEWFIK SARI 



Aj < fM < +00 such that 

Pi (S)<Di if S £ [Xi,fJ>i] 
Pi(S) > Di if S e]\i,/J,i[. 

Hence there are at most two values of S, S = Xi and S = Hi, called the break- 
even concentrations, satisfying the equation Pi(S) = Di. We adopt the convention 
/!,: = 00 if this equation has only one solution and Xi = 00 if it has no solution. 

The global analysis of this model was considered by Hsu, Hubbell and Waltman 
[5], in the Monod case [9] when the response functions are of Michaelis-Menten 
form, 

"< s » - <£h- (") 

and the yields are constant y%(S) = Yi, and Di = D for i = 1 • • -n. The authors 
showed that only the species with the lowest break-even concentration survives. 
Thus the competitive exclusion principle (CEP) holds: only one species survives, 
namely the species which makes optimal use of the resources. Hsu [4] applied a 
Lyapunov-LaSallc argument to give a simple and elegant proof of the result in [5] 
for the case of different removal rates Di. The Lyapunov function Vh used by Hsu 



is 



where 



V H = -da + ci ±-^l d £, + J2*Xi, (1-3) 



Ai 



i=2 



C = l^i-, i = l...n, and ^=^. 
Yidi-Di' 1 D l 

CEP has been proved under a variety of hypothesis [3J [51 [UJ [TS] . For a survey 
on the contribution of each paper the reader may consult the introduction of the 
papers [8j [12]. The variable yield case was considered, for n . — 1 and n = 2 
by Pilyugin and Waltman [10], with a particular interest to linear and quadratic 
yields, and by Huang, Zhu and Chang [6]. The model (jl.lj) was considered by 
Arino, Pilyugin and Wolkowicz pQ and Sari and Mazenc [12] . 

In this paper we generalize |14j by allowing variable yields and we generalize 
[T] by allowing multi species. We further extend |12| by providing less restrictive 
assumptions on the system. For biological motivations concerning the dependancc 
of the yields on the substrate, the reader is refeered to [T] [TU] and the references 
therein. 



2. Analysis of the model 

It is known (see Theorem 4.1 [1]) that the non- negative cone is invariant under 
the flow of (jl.ip and all solutions are defined and remain bounded for all t > 0. Sys- 
tem (jl.ip can have many equilibria: the washout equilibrium Eq = (S°,0, ■ ■ ■ ,0), 
which is locally exponentially stable if and only if for all i = 1 • ■ • n, S° ^ [A,,/^] 
and the equilibria E* and E** where all component of E* and E** vanish except 
for the first and the (i + l)th, which are 

S = X h Xi = x*i :=Fi(Xi), iorE* 

and 

S = fii, Xi = x** := Fi(fii), for E** 
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respectively, where 

CO _ a 

F ' (S > = D W 

The equilibrium E* lies in the non- negative cone if and only if A; < S°. If Ai < Xj 
for all i ^ j and i^' (A$) < then it is locally exponentially stable. It coalesces with 
Eq when Ai = S a . The equilibrium E** lies in the non- negative cone if and only if 
/ij < S° and is locally exponentially unstable if it exists. Its coalesces with E when 
/ii = S a . Besides these equilibria, the system (jl.ip can have a continuous set of 
non-isolated equilibria in the non-generic cases where two or more of the break-even 
concentrations are equal. In what follows we assume, that 

Ai < A 2 < • • • < A„, and Ai < S° < pb X . (2.1) 

Hence E is locally exponentially unstable and the equilibrium E± lies in the non- 
negative cone. It is locally exponentially stable if and only if 

FKM) < /i(Ai) + /((Ai)(5° - Ax) > 0. (2.2) 

We consider the global asymptotic stability of E*. 

Before presenting the results, we need the following lemma, 

Lemma 2.1. The solutions S(t), Xi(t), i = l---n of 11. 1]) with positive initial 
conditions are positive and bounded, and if Ai < S° < /Lt, for some i = 1 • • • n, then 
S(t) < S° for all sufficiently large t. 

Proof. The proof is similar to the proof of Lemma 2. 1 in [M] obtained for the model 
(|1.1|) in the case where the yields are constant. □ 

Theorem 2.2. Assume that 

(1) Ai < A 2 < ■ ■ • < \ n , and Ai < S° < (J,i. 

(2) There exist constants ai > for each i > 2 satisfying A; < S° such that 



max gi(S) < cti < min gi(S), i > 2. (2-3) 

0<S<A! \i<S<Pi 



where 



(3) The function 



MS) Pl (S)-D 1 S° -X, 

g*{b) = — — — ——^ — -. (2.4 

/i(Ai) Pi(S) - Di S° - S 



F(S) = (2.5) 



satisfies 

F(S)<F(X 1 )ifSG}0,X 1 [, F(S)>F(X 1 )ifS€]X 1 ,S°[. (2.6) 

Then the equilibrium E^ is globally asymptotically stable for system with re- 
spect to the interior of the positive cone. 



Proof. From Lemma 12.11 it follows that there is no loss of generality in restricting 
our attention to < S < S°. Consider the function V = V(S, x\, ■ • ■ , x n ) as follows 



/l(A0(S» - a) 



i=2 
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where oti are the positive constants satisfying (|2.3[) . The function V is continuously 
differentiable for < S < S° and Xi > and positive except at the point E*. The 
derivative of V along the trajectories of is given by 

v = xMS) A] F(Al ^"f (s) +Y.XMS) 

where 

h i {S) = \p i (S)-D i }[a i -g i (S)] 

and and -F^S) are given by (|2.4j) and (|2.5p respectively. First, note that, 

using (|2.1[) and (|2.6[) . the first term of the above sum is always non-positive for 
< S < S° and equals for S G]0,S*°[ if and only if S = A x or a; x = 0. If 
5 e [Ai, Ai] then < Di and p x (S) > D x so that ^ (5) < < a, for any choice 

of a, > 0. Similarly if ^ < S° and 5 € [jtt,-, 5°] then Pi(S) < Di and pi(S) > D l so 
that gi(S) <0 < on for any choice of a, > 0. On the other hand, if S G [0, Ai] then 
Pi(S) < Di and, using (j2~3|) . <?i(S) < a 4 so that /ij(S') < 0. Finally, if S 6 [Aj,p<] 
then (5) > A and gi(S) > a, so that hi(S) < 0. Thus < for every 

S e]0, provided that the numbers a l satisfy ([2~5]) . Hence V < and 7' = 
if and only if = for % = 1 • • • n or S = Ai and Xi = for i = 2 ■ • • n. We use 
LaSalle theorem (the details are as in the proof of Theorem 2.3 in [2]), the w-limit 
set of the trajectory is Ef. □ 

Condition (|2.6p is satisfied in the particular case where F'(S) > for < S < S°. 
For S = Ai we obtain the condition (|2 .2[) of local exponential stability of E'(. 

3. Applications 

In this section we show how Theorem 12.21 extends results in [TJ [TJJ E] . 

3.1. One species. In the case n = 1 the equations take the form 

S' = D(S°-S)-x 1 f 1 (S) , . 

x[ = \pi(S)- Dfa (6 - i} 

Corollary 3.1. If Xi < S° < fii and 1 — j|(f|)^gu^g] has exactly one sign change 
for S € (0, S°) then E^ is globally asymptotically stable with respect to the interior 
of the positive quadrant. 



Proof. The result follows from Theorem 12 . 2 1 since the hypothesis on the change of 
sign is is equivalent to (|2.6j) . In the case when n = 1 the condition (|2.3j) is obviously 
satisfied. □ 

Corollary 13.11 was obtained by Arino, Pilyugin and Wolkowicz (see p], Theorem 
2.11). These authors used the following Lyapunov function 



/i(Ai)(5°-a) J x . e 



Notice that the Lyapunov function (|2.7j) we use reduces to the function Vapw in 
the case when n — 1. 
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3.2. Constant yields. In the case when the yields are constant, yi(S) = Yi, the 
equations take the form 



s' = d[s° - s] - EILi 



Pi(S) 



(3.2) 

x' t = [pi(S) - D i ]x l , i = 1 ■ ■ -n, 

Corollary 3.2. Assume that 

(1) Ai < A 2 < • • • < A nj and X 1 < S° < /Lti. 

(2) There exist constants ai > for each i > 2 satisfying Xi < S° such that 

max gf L (S) < « t < min ft W£ (5), * > 2. (3.3) 

0<6<Ai \i<b<Pi 

where 

ft Di pi ( g) _ D . s o_ S - V A ) 

Then the equilibrium E\ is globally asymptotically stable for system US. 2\) with re- 
spect to the interior of the positive cone. 

Proof. In the case when the yields are constant, we have 

9i{S) = ^-gY L {S). 

Hence the conditions (|3 . 3[) imply the conditions (|2.3p with constants on replaced 
by atiYi/Yi. On the other hand F(S) = (go^gj^ • Thus, assumption (|2.6j) follows 
from hypothesis Ai < S° < The result follows from Theorem 12.21 □ 

Corollary [372] was obtained by Wolkowicz and Lu (see [14], Theorem 2.3). These 
authors used the following Lyapunov function 

VWL = / TTTT^ da + — I 7 + > — 2* 



Di(SO-a) Y 1 J X . g £^ 



=2 

Notice that the Lyapunov function (|2.7[) we use satisfies V = YiVwl m the case 
when the yields are constant. 

3.3. Variable yields. In [TJ], another Lyapunov function is proposed in the case 
when the yields are variable, leading to the following result. 

Corollary 3.3. Assume that 

(1) Ai < A 2 < • • • < A„ , and Ai < S° < 

(2) There exist constants ai > for each i > 2 smc/i i/iai 

max . gf M (S) < a,- < min gf M (S), i>2. (3.5) 

0<6<Ai \i<b<pi 



where 



9? M (S) = (3-6) 
pi (5) - A 



(3) The function F(S) = i^jg satisfies 

F(S)<F(X 1 )ifS€}0,X 1 [, F(S)>F(Xx) ifSejXuS^. (3.7) 

Then the equilibrium E\ is globally asymptotically stable for system with re- 

spect to the interior of the positive cone. 



(i 



TEWFIK SARI 



Proof. The condition (|3.7|) is the same as the condition ([2.6)1 . On the other hand 
we have 

By p.5p and ()3.7p we have (|2.3p with appropriate constants a*. The result follows 
from Theorem O □ 

Corollary 13.31 was obtained by Sari and Mazenc (see [12], Theorem 2.2). These 
authors used the following Lyapunov function 

KSM = / Tl — \ « cr + / 1 «?+> OliXi 

J Ai AW J^j £ ^ 

where Qfj , i = 2 • • ■ n are positive constants satisfying (|3.5[) . In the case when the 
response functions are of Michaelis-Menten form ([1.2)1 and the yields are constant 
the Lyapunov function Vsm reduces to the Lyapunov function Vh ([1.3|1 useb by 
Hsu [3|. Notice that the Lyapunov function ([2.7)1 we use is not proportional to the 
function Vsm- 



4. Discussion 

In this paper we considered a mathematical model ([1.1)1 of n species of microor- 
ganisms in competion in a chemostat for a single resource. The model incorporates 
both monotone and non-monotone response functions, distinct removal rates and 
variable yields. We demonstrated that the CEP holds for a large class of growth 
functions and yields. 

Even with constant yields, the problem is not yet completely solved: the CEP 
holds for a large class of growth functions [2j [3j 01 HI [14j [15| but an important open 
question remains: is the CEP true assuming only that the /, are monotone with no 
restriction on the Di ? This major open problem remains unresolved after more 
than thirty years [7]. However, in the case of constant yields numerical simulations 
of model (|1.1|) have only displayed competitive exclusion. 

In the case where the yields are constant, it is known [3] that the CEP holds 
provided that D, = D for all i, the set Q = \J ieN ]^i, fu[ is connected, and S° G Q, 
where N — {i : \i < S }. Wolkowicz and Lu [14] conjectured that this result 
can be extented to the case of different removal rates. Under condition ([2.1| . it is 
clear that the set Q is connected, and S° G Q- The condition Ai < for i ^ 1 
can be stated without loss of generality, by labelling the populations such that the 
index i = 1 corresponds to the lowest break-even concentration, but the condition 
Ai < S° < fi\ in (]2.ip cannot be stated without loss of generality. If < S°, 
it is not possible to show the CEP by the methods that we used. To the best of 
our knowledge, in the case of different removal rates and non-monotone response 
functions, the CEP has been proved only under the assumption S° < ;Ui [8[ [Ml , fT5] . 
However, Rapaport and Harmand [llj considered the case of two populations and 
proposed conditions on the growth functions such that the CEP holds under the 
condition \i\ < S°. It should be interesting to extend their methods to more general 
cases. We leave this problem for future investigations. 

Our result concern also the case of variable yields, for which it is known [H^fTO] 
that more exotic dynamical behaviours, including limit cycles, are possible. Thus 
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in the case of variable yields, it is of great importance to have criteria ensuring the 
global convergence to an equilibrium with at most one surviving species. 
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